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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE NO. 1076 

A SIMPLIFIED METHOD FOR DETERMINING PROM PLIGHT DATA 
T HE RATE OP CHANGE OP YAWING-MOMENT COEFFICIENT 

WITH SIDESLIP 

By Robert C. Bishop and Harvard Lomax 


SUMMARY 


A method is presented by which the directional stability 
derivative ®np» ra * e of change of yawing— moment coeffi- 

cient with sideslip angle, can be evaluated for a conventional 
airplane from flight records of a lateral or directional 
oscillation. For the method shown, the calculation of C n 

for a particular high— speed— flight condition reduces to the 
determination of- only the moment of inertia about the Z— axis 
and the period of a sideslipping, yawing or rolling oscilla— 
t ion. 

When applied to conventional airplanes flying at low to 
moderate lift coefficients, the assumptions Involved in this 

implified method produce negligible error. A comparison of 

ng as determined in flight and in the wind tunnel shows good 
agreement for the four conventional airplanes considered. 


INTRODUCTION 


During preliminary flight tests of experimental pursuit- 
type airplanes an accurate estimation of maximum sideslip angles 
attainable in rolling pull-outs or other dynamic conditions 
may be required before such critical tests are undertaken. 

Existing theoretical methods for calculating the maximum. __ 

sideslip angle of an airplane in dynamic flight require an 

accurate knowledge of ® n 6* Methods now used for compute— — 

tion of C a g, when wind— Tunnel data are unavailable, involve 

estimates of tail effectiveness which, because of fuselage 
interference, are subject to considerable error. Therefore a 
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method has teen developed for measuring in flight values of 
C n ^ which may he used to verify or correct the design com- 

putat ions . 

The derivation and application of equations and the 
correlation of flight results with wind-tunnel data for the 
four airplanes shown in figure 1 are presented in this report 


COEFFICIENTS AND SYMBOLS 


Coefficients and symbols defined herein are referred to 
the wind system of coordinates in which the origin is fixed 
at the center of gravity of the airplane, the Z-axis is in 
the plane of symmetry and perpendicular, to the relative air 
stream, the X— axis is in the plane of symmetry and parallel 
to the relative air stream, and the Y— axis is perpendicular 
to the Z— and X— axes. 
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wing span, feet 

wing area, square feet 

mass of airplane, slugs 

distance from center of gravity to rudder hinge line, 
feet 

acceleration due to gravity, feet per second per second 
moment of inertia about X— axis, slug-feet square 

moment of Inertia about Z— axis, slug— feet square 

air density, slugs per cubic foot 

velocity of airplane a.long flight path, feet per second 

free— stream dynamic pressure (ip? 2 ), pounds per square 
foot 

operator 

root of the stability quartic 

angle of sideslip (positive when right wing is forward) 
radians 
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3° 180 3/tt, degrees 

F period of oscillation, seconds 

p rate of roll., radians per second 

r rate of yaw, radians per second 

v component of flight velocity along the Y— axis, feet per 
second 

Y force along Y— axis, pounds 
I moment about X— axis, foot-pounds 
N moment about Z— axis, foot-pounds 
C y lateral— force coefficient (Y/qS) 

rolling-moment coefficient (L/qSb) 

C n yawing— moment coefficient (N/q.Sb) 

C Y (dCy/oP) 

P 

c n p (0C n /bP) 
c t r [bC l /b(rb/2V)] 

Cn r [bC n /b(rb/2Y)] 

C i [bC l /b(pb/2Y)] 

<V [bC n /b(pb/27)] 

Y v US/m7)C Yg 

Lg (qSb/l x )C lp 

Sp (<LSl/l Z )O n? 
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N p (o_Sb/I 2 ) (V 27) C 

ir 

N r (qSVlz)(V2V)Cn r 
L p (qS-b/lx) (V2V)0l 

Jr • • •• 

L r (qSb/l x )(b/2T)0 lr 

THEORY 

The stability quartic arising from the consideration of 
small lateral disturbances from steady horizontal flight can 
be written (reference 1) 

s 4 + Aj, s 3 + A g a 3 + Aj_ s + A 0 =0 
where 


A fe = L p (Y v + N r ) +N r Y v + N£-L r N p 

\ - “ Y v <* r \ -l r »p) - S + 1 3 C N P - <e/T)] __ 

A 0 - (e/T) (Ip K r - L r Hp) 

The real roots and the real parts of the complex roots of 
the preceding quartic determine the dapping or divergence of 
the airplane motion, and the imaginary parts determine the 
period of oscillation. Normally one of the real roots is 
small enough so that it can be approximated by neglecting 
terms in the quartic of higher order thah the first. Denoting 
this root by \ x , it follows that 

~ -u 0 /V 

By factoring this root out of the quartic the following cubic 
is obtained: 

s + (A 3 + Xj.) s + (Ajg + A 3 ) s 

+ (A x + A £ + A3 + Aj_3 ) = .0 

Eor conventional airplanes traveling at relatively high 
speeds, the values of \ x , N p , i ( Y^ — N r )S, and g/Y are 
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small enough to he neglected in comparison with the terms 
L d and Np. By these assumptions the cubic may he factored 

giving - - — — 

1 1 , 

(s-Lp) [s -- 2 (Y v + N r ) [s - 2 (Y v + N r ) + = 0 

Since Kg is positive, .yfri is the magnitude of the 
imaginary part of the complex roots and therefore the fre- 
quency of oscillation. It follows that the period of oscil- 
lation of the laterally or directionally disturbed airplane 
is given by the formula 

p « _§!!_. 

JW 


Hence 


o _ _ n fi oo 

° n * " 5.7.3 qSb 'P 2 ~ °* 688 qSb P* 

The magnitude of the error in the period introduced by 
the assumption that the period of oscillation equals 2tt 
divided by J Ng depends upon both the relative and the 
absolute magnitude of the neglejcted stability derivatives 

of the airplane. Therefore, caution should be used when 
the assumption is applied to .an unconventional design, such 
as a tailless airplane, where the relative magnitudes of 
the derivatives mey differ considerably from those of the 
airplanes of this report. 

In ord er to illustrate the accuracy of the approximation 
P = 2tt l Tg , the variations with indicated airspeed of the 
period as computed by three methods are shown in figure 2 
for a representative modern airplane. The three methods used 
are: (1) the method of reference 1 in which (within the 

limitations of the initial assumptions) the theoretically 
true period is given as 2tt divided by the magnitude of the 
complex part of an imaginary root of the stability quartic, 

(2) the method of the present report in which the period 
is equal td 2tt / vlTp, and (3) the method of reference 2 in 
which an alternative simplification and factorization of the 
quartic gives a period equal to . ... 


2tt 



Lb (n p -g/v) 

-Nr — Lip 


zleJb 
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The comparisons of figure 2 indicate that, for the airplane 
considered, the approximation of this report (method (2)) 
gives results that agree within 5 percent with those obtained 
using the theoretically true method (method (1)), and are 
closer than those of method (3). At low speeds (and high 

lift coefficients), where the neglected stability deriva- 

tives assume larger values relative to Cng, the approxi- 
mation presented herein tends to become less reliable and 
others should be used. 


PLIGHT PROCEDURE 


The period of motion necessary for calculation of C n p 
may be obtained by placing the airplane in a steady sideslip 
while trimmed for straight and zer o— s ide s lip flight, then 
abruptly returning the control to trim and recording the 
resulting motion. A typical time history of the yawing and 
rolling oscillation that follows such a disturbance is shown 
in figure 3. Prom records similar to those shown in figure 3 
the period of the oscillations may be determined and inserted 
in the equation given for Onp . However, it is net essential 
in the evaluation of C n g to obtain time histories of the 

airplane motions; any accurate means of measuring the period_ 
alone would be sufficient. 

The value of C n g obtained by this method actually 
represents an average value for the range of sideslip angles 
covered. If it is desired to establish a curve of C n g 

against 6, runs can be made with varying amounts of applied 
disturbance or with different initial trim positions. 

The procedure moy be used for either a control— fixed or 
a control— free configuration. In the interpretation of 

control— free results consideration should be given to the 

control— system friction, which is sometimes sufficient to 
hold the rudder fixed in a rudder— free maneuver. 


PRECISION 

The accuracy with which C n can bo evaluated from 
flight tests depends largely cn the precision with which t he 
period can be measured and the moment of inertia Ig can be 
estimated. If e represents the error in the flight— measured 
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period, the percent error in Cng resulting from the error _ 
e is given by 

Percent error in On- = ±100 * 2e/P 

P 

It is estimated on the basis of available flight data that 
the value of £ is of the order of 0,10 second. In order 
to reduce the percentage error from this source to less than 
6 percent in any given reading, a period greater than about 
3.5 seconds would be required. Since the period is inversely 
proportional to the speed, the experimental error is reduced 
for lower speeds. (The error likely to be incurred in the 
computation of I^ is greatly dependent upon the method used 
for its determination and with care may be held within 1 
percent, 

The gain in experimental accuracy at lower speeds tends 
to be offset by the decreased accuracy in the theory due to 
the approximations. It is expected that the greatest over— all 
accuracy will be obtained between 200 and 300 miles per hour. 


COMPARISON OP PLIGHT DATA WITH WIND-TUNNEL DATA 

To check the validity of the simplified method developed 
in this report a comparison has been made of the values of 
C n g obtained from wind-tunnel and flight tests. Pour air- 
planes were chosen (fig. 1) which exhibited a wide range of 
directional stability and for which there existed the 
necessary wind-tunnel and flight data. 

The value of C no was obtained from flight data using 
/ 0 r 

C np0 = 0.6.8 8 “g^pps — 

The moments of inertia were obtained from reports of the 
Langley Spin Tunnel, and the values of P were determined 
from directional oscillations where the amplitude of the side- 
slip oscillations was ± 5°. The flight and wind-tunnel results 
in the form of C n against are compared in figure 4. 

Since no values of intercept can be assigned to this curve, 
the flight curves were arbitrarily placed so that zero inter- 
cepts coincided with those of the wind-tunnel data. Good 
agreement between wind-tunnel and flight— test results for all_. 
four airplanes indicates that results at least as accurate as 
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those obtainable in a wind tunnel may be anticipated if the 
approximate method of this report is used. 


CONCLUDING REMARKS 


The simplified method for the determination of 
by flight te.st requires the measurement of only the period 
of oscillation of the airplane when disturbed laterally or 
directionally. lor conventional airplanes flying at low to 
moderate lift coefficients, the assumptions involved in this 
simplified method produce negligible error; hence tlpe method 
gives results which are more accurate than those obtained 
using the approximate method of reference 2 and which are in 
good agreement with wind— tunnel results. For unconventional 
airplanes (such as tailless) or for airplanes flying at high 
lift coefficients the error nay become appreciable and the 
method should be used with caution . 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., March 1946. 
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Figure 2.- Comparison of the true and approximate methods for 
determining the theoretical period of a representa- 
tive airplane. Stability derivatives were computed for a lift 
coefficient of 0.132. 
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Fig. 3 



Figure 3.- Time history of a typical rudder-fixed lateral 
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Fig. 4 
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Figure 4.- Comparison of flight and wind-tunnel values 
of rate of change of yawing-moment coeffic- 
ient with sideslip angle for four airplanes. 




